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We report first-principles density-functional theory studies of native point defects and defect 
complexes in olivine-type LiFePCU, a promising candidate for rechargeable Li-ion battery electrodes. 
The defects are characterized by their formation energies which are calculated within the GGA+?7 
framework. We find that native point defects are charged, and each defect is stable in one charge 
state only. Removing electrons from the stable defects always generates defect complexes containing 
small hole polarons. Defect formation energies, hence concentrations, and defect energy landscapes 
are all sensitive to the choice of atomic chemical potentials which represent experimental conditions. 
One can, therefore, suppress or enhance certain native defects in LiFePCU via tuning the synthesis 
conditions. Based on our results, we provide insights on how to obtain samples in experiments with 
tailored defect concentrations for targeted applications. We also discuss the mechanisms for ionic and 
electronic conduction in LiFePCU and suggest strategies for enhancing the electrical conductivity. 
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I. INTRODUCTION 



Olivine-type LiFeP04 is a promising candidate for 
rechargeable Li-ion battery electrodes^ The material is 
known for its structural and chemical stabilities, high in- 
tercalation voltage (^3.5 V relative to lithium metal), 
high theoretical discharge capacity (170 mAh/g), envi- 
ronmental friendliness, and potentially low costs^^ The 
major drawback of LiFeP04 is poor ionic and electronic 
conduction (with an electrical conductivity of about 10~ 9 
S/cm at 298 K)i that limits its applicability to devices. 
While the conduction can be improved by, e.g., mak- 
ing LiFeP04 nanoparticles and coating with conductive 
carbon^i^ the high processing cost associated with the 
manufacturing of carbon-coated LiFeP04 nanoparticles 
may make it less competitive than other materials. An- 
other approach is to dope LiFeP04 with aliovalent im- 
purities (Mg, Ti, Zr, Nb) which was reported to have en- 
hanced the conductivity by eight orders of magnitude. 7 
The role of these dopants in the conductivity enhance- 
ment, however, remains controversial^ 9 - A better under- 
standing of aliovalent doping, and also better solutions 
for improving the performance, first requires a deeper 
understanding of the fundamental properties, especially 
those associated with native defects, which is currently 
not available. First-principles density-functional theory 
(DFT) studies of native point defects and defect com- 
plexes in LiFeP04 can help address these issues. 

It is now generally accepted that LiFeP04 is an insu- 
lating, large band-gap material in which electronic con- 
duction proceeds via hopping of small hole polaronsJ^— 
These polarons may be coupled to other defects such as 
lithium vacancies j 11 ' 12 Iron antisites (Feu) have also been 
reported to be present in LiFeP04 samples*^ - — This na- 
tive defect is believed to be responsible for the loss of elec- 
trochemical activity in LiFeP04 due to the blockage of 
lithium channels caused by its low mobility] 18 ' 19 Clearly, 



native defects have strong effects on the material's perfor- 
mance. Experimental reports on the defects have, how- 
ever, painted different pictures. Some authors reported 
evidence of some iron and lithium atoms exchanging sites 
and forming the antisite pair FeLi-Lipe) 17 ' 19 while oth- 
ers determined that Feu is formed in association with 
lithium vacancies (VlO' 15 ' 18 These conflicting reports 
suggest that the results may be sensitive to the actual 
synthesis conditions, and indicate that a better under- 
standing of the formation of native defects in LiFeP04 
is needed in order to produce samples with controlled 
defect concentrations. 

Computational studies of native defects in LiFeP04 
and related compounds have been reported by several 
research groups i 11 ^ " — Notably, Maxisch et al. studied 
the migration of small hole polarons in LiFeP04 using 
first-principles calculations where the polarons were cre- 
ated both in the absence and in the presence of lithium 
vacancies^ The first systematic study of native de- 
fects in LiFeP04 was, however, carried out by Islam 
et al. using interatomic-potential simulations where they 
found the antisite pair FeLi-Lipe to be energetically most 
favorable i 21 i 22 Based on results of first-principles calcu- 
lations, Malik et al. recently came to a similar conclusion 
about the antisite paired Although these studies have 
provided valuable information on the native defects in 
LiFeP04, they have three major limitations. First, stud- 
ies that make use of interatomic potentials may not well 
describe all the defects in LiFeP04. Second, these stud- 
ies seem to have focused on neutral defect complexes and 
did not explicitly report the structure and energetics of 
native point defects as individuals. Third, and most im- 
portantly, none of these previous studies have thoroughly 
investigated the dependence of defect formation energies 
and hence defect concentrations on the atomic chemical 
potentials which represent experimental conditions dur- 
ing synthesis. 
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We herein report our first-principles studies of the 
structure, energetics, and migration of native point de- 
fects and defect complexes in LiFePC>4. We find that 
defect formation is sensitive to the synthesis conditions. 
Native defects can occur in the material with high con- 
centrations and therefore are expected to have impor- 
tant implications for ionic and electronic conduction. We 
will show how conflicting experimental data on the native 
defects can be reconciled under our results and provide 
general guidelines for producing samples with tailored 
defect concentrations. Comparison with previous theo- 
retical works will be made where appropriate. In the 
following, we provide technical details of the calculations 
and present the theoretical approach. Next, we discuss 
the structural and electronic properties of LiFeP04 which 
form the basis for our discussion of the formation of na- 
tive defects in the material. We then present results of 
the first-principles calculations for native point defects 
and defect complexes, focusing on their formation ener- 
gies and migration barriers, and discuss the dependence 
of defect formation energies on the atomic chemical po- 
tentials. Based on our results, we discuss the implications 
of native defects on ionic and electronic conduction, and 
suggest strategies for enhancing the electrical conductiv- 
ity. Finally, we end this Article with some important 
conclusions. 



II. METHODOLOGY 

Computational Details. Our calculations were 
based on density- functional theory within the GGA+U 
framework,— ~— which is an extension of the generalized- 
gradient approximation (GGA),— and the projector 
augmented wave method ) 29 ' 30 as implemented in the 
VASP code^i^ In this work, we used £7=5.30 eV and 
J=1.00 eV for iron in all the calculations (except oth- 
erwise noted), i.e., the effective interaction parameter 
[7— , 7=4.30 eV (hereafter U—J will be referred to as U 
for simplicity). This value of U is the averaged value 
based on those Zhou et al. calculated self-consistently for 
iron in LiFeP0 4 (i.e., Fe 2+ : [7=3.71 eV) and in FeP0 4 
(i.e., Fe 3+ : [7=5.90 eV), which has been shown to cor- 
rectly reproduce the experimental intercalation poten- 
tial of LiFeP04.— It is known that the results obtained 
within GGA+C7 depend on the value of U. However, we 
have checked the U dependence in our calculations and 
find that the physics of what we are presenting is insen- 
sitive to the U value for 3.71 eV < U < 5.90 eV. 

Calculations for bulk olivine-type LiFeP04 (or- 
thorhombic Pnma; 28 atoms/unit cell) were performed 
using a 4x7x9 Monkhorst-Pack k-point mesh.— For de- 
fect calculations, we used a (1x2x2) supercell, which cor- 
responds to 112 atoms/cell, and a 2x2x2 k-point mesh. 
The plane-wave basis-set cutoff was set to 400 eV. Con- 
vergence with respect to self-consistent iterations was 
assumed when the total energy difference between cy- 
cles was less than 10 -4 eV and the residual forces were 



less than 0.01 eV/A. In the defect calculations, the lat- 
tice parameters were fixed to the calculated bulk val- 
ues, but all the internal coordinates were fully relaxed. 
The migration of selected defects in LiFeP04 was stud- 
ied using the climbing-image nudged elastic-band method 
(NEB). 36 All calculations were performed with spin po- 
larization and, unless otherwise noted, the antiferromag- 
netic spin configuration of LiFePCU was usedi^l 

Defect Formation Energies. Throughout this Ar- 
ticle, we employ defect formation energies to character- 
ize different native point defects and defect complexes in 
LiFeP04. The formation energy of a defect is a crucial 
factor in determining its concentration. In thermal equi- 
librium, the concentration of the defect X at temperature 
T can be obtained via the relatio n 38 ' 39 

c(X) = A sitcs A config exp[-£'(X)/fc B T], (1) 

where N s i tcs is the number of high-symmetry sites in the 
lattice per unit volume on which the defect can be in- 
corporated, and iVconfig is the number of equivalent con- 
figurations (per site). Note that the energy in Eq. ([l} 
is, in principle, a free energy; however, the entropy and 
volume terms are often neglected because they are neg- 
ligible at relevant experimental conditions^ 9 - It emerges 
from Eq. (JXJ) that defects with low formation energies will 
easily form and occur in high concentrations. 

The formation energy of a defect X in charge state q 
is defined as^ 8 - 

£^(X«) = £ t ot(X 9 ) - £ t ot(bulk) -J2 n ^ 

i 

+q(E v + AV + e F ), (2) 

where E tot (X. q ) and -Etot(bulk) are, respectively, the to- 
tal energies of a supercell containing the defect X and of 
a supercell of the perfect bulk material; fii is the atomic 
chemical potential of species i (and is referenced to the 
standard state), and Hi denotes the number of atoms 
of species i that have been added (rii>0) or removed 
(rii<0) to form the defect. e F is the electron chemical 
potential, i.e., the Fermi level, referenced to the valence- 
band maximum in the bulk (E v ). AV is the "potential 
alignment" term, i.e., the shift in the band positions due 
to the presence of the charged defect and the neutraliz- 
ing background, obtained by aligning the average elec- 
trostatic potential in regions far away from the defect to 
the bulk valued 8 - Note that we denote defect X in charge 
state q as X 9 . For example, Fe^; indicates that defect Feu 
occurs with charge q—+l, which is equivalent to Fe*j in 
the Kroger- Vink notation. For a brief discussion on the 
use of notations, see, e.g., Ref4°-. 

Chemical Potentials. The atomic chemical poten- 
tials fii are variables and can be chosen to represent ex- 
perimental conditions. \ii can, in principle, be related to 
temperatures and pressures via standard thermodynamic 
expressions. The chemical potential for O2 in oxygen gas, 
for example, is given by 41 - 

Po 2 ( T ,P) = Ho 2 {T,Po) + fcTln— , (3) 

Po 
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where p and p are, respectively, the partial pressure and 
reference partial pressure of oxygen; k is Boltzmann's 
constant. This expression allows us to calculate fio 2 (T, p) 
if we know the temperature dependence of no 2 (T,p ) at 
a particular pressure p . In this work, we choose the 
reference state of ^o 2 (T, p) to be the total energy of an 
isolated 2 molecule (Eg*)& 

The value of fXi is subject to various thermodynamic 
limits. For LiFePC>4, the stability condition requires that 

A*Li + A*Fo + MP + 2/zo 2 = Aff / (LiFeP0 4 ), (4) 

where AH? is the formation enthalpy. This condition 
places a lower bound on the value of fi{. Additionally, 
one needs to avoid precipitating bulk Li, Fe, and P, or 
forming O2 gas. This sets an upper bound on the chem- 
ical potentials: /ii<0<2S. There are, however, further con- 
straints imposed by other competing Li-Fe-P-02 phases 
which usually place stronger bounds on /Xj. For example, 
in order to avoid the formation of Li3P04, 

+ W + 2^o 2 < Ai^(Li 3 P0 4 ). (5) 

After taking into account the constraints imposed by 
all possible competing phases, one can define the chem- 
ical potential range of Li, Fe, and O2 that stabilizes 
LiFeP04 which is, in fact, bound in a polyhedron in 
the three-dimensional (nu, HFc, Ho 2 ) space. For a given 
point in the polyhedron, one can determine the remain- 
ing variable /Ltp via Eq. (|4|). In this work, the formation 
enthalpies of all different Li-Fe-P-02 phases are taken 
from Ong et al4^ who have computed the energies using 
a methodology similar to ours. For example, the calcu- 
lated formation enthalpy of LiFePC>4 at T=0 K (with re- 
spect to its constituents) is —18.853 eV per formula unit 
(f.u.)^i almost identical to that (of -18.882 eV/f.u.) ob- 
tained in our calculations. Ong et al. have also calculated 
the phase diagrams of the quaternary Li-Fe-P-02 system 
at K that involve all possible phases between Li, Fe, 
P, and O2. These phase diagrams show LiFeP04 is sta- 
ble over a range of the oxygen chemical potential values, 
from —11.52 eV (where the first Fe 2+ -containing phase 
appears) to —8.25 eV (the last of the Fe 2+ -containing 
phosphates being reduced) 4^ This corresponds to /io 2 
ranging from —3.03 to —8.25 eV with respect to our cho- 
sen reference (Eq*). 

Figure [1] shows the slice of the (nu, MFc, Mo 2 ) polyhe- 
dron in the iio 2 —— 4.59 eV plane, constructed with the 
calculated formation enthalpies (taken from Ref.— ) for 
different Li-Fe-P-02 phases. The shaded area (marked 
by Points A, B, C, D, and E) shows the range of /UTj and 
/XFe values where LiFeP04 is stable. Point A, for exam- 
ple, corresponds to equilibrium of LiFeP04 with Fe2 03 
and Fe3(P04)2- At this point in the chemical-potential 
diagram, the system is close to forming Fe-containing 
secondary phases (i.e., Fe203 and Fe 3 (P04)2) and far 
from forming Li-containing secondary phases. This can 
be considered as representing a "Li-deficient" environ- 
ment. Similarly, Point D can be considered as represent- 
ing a "Li-excess" environment, where the system is close 
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FIG. 1: Chemical-potential diagram for LiFePCU at 
fJ,o 2 =— 4.59 eV. The fio 2 axis extends out of the page. Only 
phases that can be in equilibrium with LiFeP04 are included 
and the lines delineating these phases define the stability re- 
gion of LiFeP04, here shown as a shaded polygon. 

to forming Li-containing secondary phases (i.e., Li4P207 
and Li 3 P04). Note that "Li-deficient" and "Li-excess" 
environments in this sense do not necessarily mean that 
/iLi in the latter is higher than in the former, as seen in 
Fig. [U Reasonable choices of the atomic chemical poten- 
tials should be those that ensure the stability of the host 
compound. In the next sections we will present our cal- 
culated formation energies for various native defects in 
LiFeP04 and discuss how these defects are formed under 
different experimental conditions. 

Defect Complexes. Native point defects in LiFeP04 
may not stay isolated but could instead agglomerate and 
form defect complexes. For a complex XY consisting of 
X and Y, its binding energy Et> can be calculated us- 
ing the formation energy of the complex and those of its 
constituents^ 

E b = E f (X) + E s (Y) - E f (XY) , (6) 

where the relation is defined such that a positive binding 
energy corresponds to a stable, bound defect complex. 
Having a positive binding energy, however, does not mean 
that the complex will readily form. For example, under 
thermal equilibrium, the binding energy Eb needs to be 
greater than the larger of E? (X.) and E^(Y) in order 
for the complex to have higher concentration than its 
constituents^ For further discussions on the formation 
of defect complexes, see, e.g., Ref<2£. 

III. BULK PROPERTIES 

Before presenting our results for native defects in 
LiFeP04, let us discuss some basic properties of the pris- 
tine compound. Olivine-type LiFeP04 was reported to 
crystallize in the orthorhombic space group Prima with 
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FIG. 2: Electronic density of states (DOS) of LiFeP0 4 in (a) 
antiferromagnetic (AFM) and (b) ferromagnetic (FM) spin 
configurations. The zero of the energy is set to the highest 
occupied state. 



a=10.3377(5), 6=6.0112(2), and c=4.6950(2) A.- 37 The 
compound can be regarded as an ordered arrangement 
of Li + , Fc 2+ , and (P04) 3 ~ units. Li + forms Li channels 
along the 6-axis whereas Fe 2+ stays at the center of a 
slightly distorted FeC>6 octahedron (interwoven with PO4 
tetrahedra). This simple bonding picture will be very 
useful when interpreting the structure and energetics of 
native defects in LiFePCj. The calculated lattice param- 
eters are a=10.461, 6=6.061, and c=4.752 A, in satisfac- 
tory agreement with the experimental values. The cal- 
culated values are slightly larger than the experimental 
ones as expected since it is well known that GGA tends 
to overestimate the lattice parameters. The calculated 
magnetic moment for iron (Fe 2+ ) is 3.76 //b, comparable 
to the experimental value of 4.19(5) /ib at 2 Ki 3 - 7 - 

Figure [2] shows the total electronic density of states 
of LiFePC>4 in antiferromagnetic (AFM) and ferromag- 
netic (FM) spin configurations. An analysis of the wave- 
functions shows that, in both configurations, the valence- 
band maximum (VBM) and conduction-band minimum 
(CBM) are Fe 3d states. Between the highly localized 
Fe d states just below the Fermi level (at eV) and the 
lower valence band (which consists predominantly of O 
2p and Fe 3d states) there is an energy gap of about 0.40 
eV (AFM). The Li 2s state is high up in the conduction 
band, suggesting that Li donates its electron to the lat- 
tice and becomes Li + . There is strong mixing between P 
3p and O 2p states, indicating covalent bonding within 
the (P04) 3 ~ unit. The calculated band gap is 3.62 and 
3.58 eV for AFM and FM spin configurations, respec- 
tively, in agreement with previously reported value (of 
3.7 eV)iiS Experimentally, LiFeP04 has been reported 
to have a band gap of about 3.8—4.0 eV, obtained from 
diffuse reflectance measurements! 10 i 13 The compound is 
therefore an insulating, large band-gap material. 

In the GGA+[7 framework, the electronic structure 
can depend on the U value. Indeed, we find that the 
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FIG. 3: Calculated formation energies of native point de- 
fects and defect complexes in LiFeP04, plotted as a func- 
tion of Fermi level with respect to the VBM. The energies 
are obtained at Point A in the chemical-potential diagram 
for /xo 2 = — 4.59 eV (cf. Fig. [T]), representing equilibrium with 
Fe 2 3 and Fe 3 (P0 4 ) 2 . 



calculated band gap of LiFeP04 is 3.20 and 4.00 eV in 
the AFM spin configuration for £7=3.71 and 5.90 eV, 
respectively, compared to 3.62 eV obtained in calcula- 
tions using [7=4.30 eV mentioned earlier. The energy 
gap between the highest valence band (Fe 3d states) and 
the lower valence band (predominantly O 2p and Fe 3d 
states) is also larger for smaller [7 value: 0.58 and 0.20 
eV for [7=3.71 and 5.90 eV, respectively. However, our 
GGA+U calculations show that the electronic structure 
near the band gap region is not sensitive to the choice 
of [7 value, for [7 lying within the range from 3.71 to 
5.90 eV. As we illustrate in the next section, knowing 
the structural and electronic properties, especially the 
nature of the electronic states near the VBM and CBM, 
is essential in understanding the formation of native de- 
fects in LiFeP04. 



IV. FORMATION OF NATIVE DEFECTS 

In insulating, large band-gap materials such as 
LiFeP04, native point defects are expected to exist in 
charged states other than neutral, and charge neutrality 
requires that defects with opposite charge states coexist 
in equal concentrations £2r— We therefore investigated 
various native defects in LiFeP04 in all possible charge 
states. These defects include hole polarons (hereafter 
denoted as p + ), lithium vacancies (VlO and interstitials 
(Lij), iron antisites (Feu), lithium antisites (Lip c ), iron 
vacancies (Vp e ), and PO4 vacancies (Vpo 4 )- We also con- 
sidered defect complexes that consist of certain point de- 
fects such as FeLi-VLi (a complex of Feu and Vu), Feu- 
Lipe (a complex of Feu arid Lip e ), and 2FeLi-Vp e (a com- 
plex of two FeLi and one Vpe). 
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Figure [3] shows the calculated formation energies of 
relevant native point defects and defect complexes in 
LiFePCU for a representative oxygen chemical potential 
value, no 2 — ~ 4-59 eV, and nu=— 2-85 eV, ^f c =— 2.18 
eV, and /zp=— 4.64 eV. This set of atomic chemical po- 
tentials corresponds to Point A in Fig. [I] representing the 
limiting case (Li-deficient) where Fe2C>3, Fe3(PC>4)2, and 
LiFeP04 are in equilibrium. The slope in the formation 
energy plots indicates the charge state. Positive slope 
indicates that the defect is positively charged, negative 
slope indicates the defect is negatively charged. With the 
chosen set of atomic chemical potentials, the positively 
charged iron antisitc FejK and negatively charged lithium 
vacancy (VjJ) have the lowest formation energies among 
the charged point defects for a wide range of Fermi-level 
values. While there are different charged point defects co- 
existing in the system with different concentrations, the 
ones with the lowest formation energies have the high- 
est concentrations and are dominant.— ~— Figure |3] in- 
dicates that, in the absence of electrically active impu- 
rities that can affect the Fermi-level position, or when 
such impurities occur in much lower concentrations than 
charged native defects, the Fermi level will be pinned at 
£i?=1.06 eV, where the formation energies and hence, ap- 
proximately, the concentrations of FeJ; and V L ~ are equal. 
Also, charged native defects have positive formation en- 
ergies only near e^=1.06 eV. Therefore, any attempt to 
deliberately shift the Fermi level far away from this po- 
sition and closer to the VBM or CBM, e.g., via doping 
with acceptors or donors, will result in positively or neg- 
atively charged native defects having negative formation 
energies, i.e., the native defects will form spontaneously 
and counteract the effects of doping.— i 39 i 46 i 47 This indi- 
cates that LiFeP04 cannot be doped p-type or n-type. 
In the following, we analyze in detail the structure and 
energetics of the native defects. The dependence of de- 
fect formation energies on the choice of atomic chemical 
potentials will be discussed in the next section. 

Small Hole Polarons. The creation of a free posi- 
tively charged (hole) polaron p + (i.e., p + in the absence 
of other defects or extrinsic impurities) involves remov- 
ing one electron from the LiFeP04 supercell (hereafter 
referred to as "the system"). This results in the forma- 
tion of a Fe 3+ site in the system. The calculated mag- 
netic moment at this (Fe 3+ ) site is 4.28 /ib, compared to 
3.76 /iB at other iron (Fe 2+ ) sites. The local geometry 
near the Fe 3+ site is slightly distorted with the neigh- 
boring O atoms moving toward Fe 3+ ; the average Fe-0 
bond length is 2.07 A, compared to 2.18 A of the other 
Fe-0 bonds. Note that in pristine FePC>4, the delithi- 
ated phase of LiFeP04, the calculated magnetic moment 
is 4.29 hb at the iron (Fe 3+ ) sites, and the calculated 
average Fe-0 bond length is 2.06 A. This indicates that 
a hole (created by removing an electron from the sys- 
tem) has been successfully stabilized at one of the iron 
sites and the lattice geometry is locally distorted, giving 
rise to a hole polaron in LiFeP04. Since the local dis- 
tortion is found to be mostly limited to the neighboring 



O atoms of the Fe 3+ site, this hole polaron is considered 
as small polaron where the hole is "self-trapped" in its 
own potentia h 48 ' 49 The formation of free hole polarons in 
LiFePC>4 is necessarily related to the rather strong inter- 
action between Fe 3c? and O p states, and the fact that 
the VBM consists predominantly of the highly localized 
d states. 

We have investigated the migration path of p + and esti- 
mated the energy barrier using the NEB method. 36 The 
migration of p + involves an electron and its associated 
lattice distortion being transferred from a Fe 2+ site to a 
neighboring Fe 3+ site. Since spin conservation is required 
in this process, we carried out our calculations not using 
the ground-state AFM structure of LiFeP0 4 but the FM 
one where all the spins are aligned in the same direction. 
We calculated the migration path by sampling the atomic 
positions between ground-state configurations. For those 
configurations other than ground-state ones, the atomic 
positions were kept fixed and only electron density was 
relaxed self-consistently, similar to the method presented 
in Ref^ 1 -. The migration barrier is the energy difference 
between the highest-energy configuration and the ground 
state. We find that the migration barrier of p + is 0.25 eV 
between the two nearest Fe sites approximately in the b-c 
plane, which is comparable to that (0.22 eV) reported by 
in Refcii. 

Vacancies and Interstitials. Negatively charged 
lithium vacancies (V L ~j") are created by removing a Li + 
ion from the system. Since, in LiFePC>4, Li donates one 
electron to the lattice one expects that the removal of Li + 
causes only a small disturbance in the system. Indeed we 
see that lattice relaxations around the void formed by the 
removed Li + are negligible. The energy needed to form 
should also be small, consistent with our results in 
Fig. Vj^iJ on the other hand, is created by removing 
a Li atom (i.e., Li + and an electron) from the system. 
This leads to the formation of a void (at the site of the 
removed Li + ) and an Fe 3+ (formed by the removed elec- 
tron) at the neighboring Fe site. Similar to the free hole 
polaron, the neighboring O atoms of the Fe 3+ site in 
also move toward Fe 3+ , with the average Fe-0 distance 
being 2.07 A. The calculated magnetic moment is 4.29 fj,B 
at the Fe 3+ site, equal to that at the Fe 3+ site in the case 
of a free polaron. V^, therefore, should be regarded as a 
complex of V L T and p + , with the two defects being 3.26 A 
apart. Figure BJa) shows the structure of Vj-J. The bind- 
ing energy of Vfj is 0.34 eV (with respect to V L 7 and p + ). 
Note that this value is 0.42 eV in our calculations using 
(1x3x3) supercells which have 252 atoms/cell. Our esti- 
mated binding energy is thus comparable to that of 0.39 
and about 0.50 eV reported by Fisher et ali^ and Max- 
isch et al.j 1 - 1 - respectively. For lithium interstitials, the 
stable defect is Li^ 1 ", created by adding Li + into the sys- 
tem. Other charge states of Vu and LL, are not included 
in Fig.[3]because they either have too high energies to be 
relevant or are unstable. 

The migration path of VJT is calculated by moving a 
Li + unit from a nearby lattice site into the vacancy. The 
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FIG. 4: Defects in LiFePCU: (a) V® { can be regarded as a complex of VjJ (represented by the empty sphere) and hole polaron 
p + (i.e., Fe 3+ ; decorated with the square of the wavefunctions of the lowest unoccupied state in the electronic structure of 
LiFePCU in the presence of V°i); (b) Fejh-VJJ, a complex of FeJ; and V£l; and (c) Fe^-Lip,,, a complex of FeJ; and Lip c . Large 
(gray) spheres are Li, medium (blue) spheres Fe, small (yellow) spheres P, and smaller (red) spheres O. 



energy barrier for is estimated to be 0.32 eV along 
the 6-axis and 2.27 eV along the c-axis. This suggests 
that, in the absence of other native defects and extrinsic 
impurities, lithium diffusion in LiFcPC>4 is highly one- 
dimensional along the Li channels (6-axis) because the 
energy barrier to cross between the channels is too high. 
The migration path of is, however, not a straight 
line but a curved path along the 6-axis. Our results are 
thus in general agreement with previously reported theo- 
retical studies 2 ^— and experimental observation^ The 
estimated energy barriers for the migration of along 
the b and c axes are lower than those (0.55 and 2.89 eV, 
respectively) reported by Islam et al—i obtained from 
calculations using inter-atomic potentials, but closer to 
those (0.27 and about 2.50 eV) reported by Morgan et 
alj2£ obtained in GGA calculations with smaller super- 
cells. For Vl*;, a complex of p + and Vj^, one can estimate 
the lower bound of the migration barrier by taking the 
higher of the migration energies of the constituents, 45 
which is 0.32 eV (along the 6-axis), the value for V{^. 

Other possible vacancies in LiFeP04 are those associ- 
ated with Fe 2+ and (P0 4 ) 3 ~ units. The creation of V§~ 
corresponds to removing Fe 2+ from the system. We find 
that this negatively charged defect causes significant re- 
laxations in the lattice geometry. The neighboring Li + 
ions move toward the defect, resulting in the Li channels 
being bent near V F ~ where Li + ions are displaced up 
to 0.27 A from their original positions. can be con- 
sidered as a complex of and p + with the distance 
between the two defects being 3.81 A. V^~, on the other 
hand, corresponds to removing Fe 2+ but leaving an elec- 
tron in the system. This defect can be regarded as a com- 
plex of Vp~ and a negatively charged (electron) polaron 
(hereafter denoted as p~ ). At the Fe site where the elec- 
tron polaron resides, which is 7.68 A from the vacancy, 
the calculated magnetic moment is 2.86 /ib; the average 
Fe-0 distance is 2.30 A, which is larger than that associ- 
ated with other Fe sites (2.18 A). Finally, Vpo 4 is stable 
as VpQ^ as expected. This positively charged defect cor- 
responds to removing the whole (PO4) 3- unit from the 



system. With the chosen set of atomic chemical poten- 
tials, Vp~ and VpQ^ have very high formation energies 
(2.33 eV and 3.56 eV, respectively, at ei?=1.06 eV) and 
are therefore not included in Fig. [3] 

Antisite Defects. Lithium antisites Lip e are created 
by replacing Fe at an Fe site with Li. Lip e can be consid- 
ered as replacing Fe 2+ with Li+. Due to the Coulombic 
interaction, the two nearest Li + ion neighbors of Lip e are 
pulled closer to the negatively charged defect with the 
distance being 3.25 A, compared to 3.32 A of the equiv- 
alent bond in pristine LiFeP04. Li Fc , on the other hand, 
can be regarded as a complex of Lip e and p + with the dis- 
tance between the two defects being 3.98 A. The binding 
energy of Li Fc (with respect to Lip c and p + ) is 0.30 eV. 
Similarly, one can replace Li at an Li site with Fe, which 
creates an iron antisite Feu- Fe^ corresponds to replac- 
ing Li + with Fe 2+ , whereas Fe 2 / can be regarded as a 
complex of Fe^; and p + . For Fe^; , which corresponds to 
replacing one Li + with Fe 2+ and adding an extra elec- 
tron to the system, the extra electron is stabilized at the 
substituting Fe atom, where the calculated magnetic mo- 
ment is 2.95 fiB- One might also regard Fe^ as a complex 
of Fe^i and p~ , but in this case the two defects stay at 
the same lattice site. With the chosen set of chemical 
potentials, Fe"; has a very high formation energy (2.04 
eV) and is therefore not included in Fig. [3] Again, other 
native defects that are not included here are unstable or 
have too high formation energies to be relevant. 

Defect Complexes. From the above analyses, it is 
clear that defects such as p + (p~), V L 7, V F ~ , FejK, Lip e , 
and Vp 3 + 4 can be considered as elementary native defects 
in LiFeP04, i.e., the structure and energetics of other 
native defects can be interpreted in terms of these basic 
building blocks. This is similar to what has been ob- 
served in complex hydrides. — These elementary defects 
(except the free polarons) are, in fact, point defects that 
are formed by adding and/or removing only Li+, Fe 2+ , 
and (P04) 3 ~ units. They have low formation energies 
(c/. Fig. [3]) because the addition/removal of these units 
causes the least disturbance to the system, which is con- 
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sistent with the simple bonding picture for LiFeP04 pre- 
sented in the previous section. The identification of the 
elementary native defects, therefore, not only helps us 
gain a deeper understanding of the structure and ener- 
getics of the defects in LiFcP04 but also has important 
implications. For example, one should treat the migra- 
tion of defects such as V^i as that of a V-T and p + complex 
with a finite binding energy, rather than as a single point 
defect. 

In addition to the defect complexes that involve p + and 
p~ such as Vl ;, Vp~, Vy~, Li Fo , Fe^ ; , and Fej"+ described 
above, we also considered those consisting of V{^, Fe^, 
Li Fc , and V$~ such as Fe^-Lip e , Fejh-V L T, and 2Fe£ r 
V3~. Figure HJb) shows the structure of FeT^-V^. The 
distance between FeT^ and is 2.96 A (along the b- 
axis), compared to 3.03 A between the two Li sites in 
pristine LiFeP04. We find that this complex has a for- 
mation energy of 0.36—0.56 eV for reasonable choices of 
atomic chemical potentials, and a binding energy of 0.49 
eV. With such a relatively high binding energy, even 
higher than the formation energy of isolated Fe^ and 
Vl7 (0.42 eV at e F =1.06 eV, cf. Fig. Fe^-V^ is ex- 
pected to occur with a concentration larger than either 
of its constituents under thermal equilibrium conditions 
during synthesis*^ In FeT^-V^, the energy barrier for mi- 
grating Fe^i to Vj^ is about 0.74 eV, comparable to that 
(0.70 eV) reported by Fisher et al32. This value is twice 
as high as the migration barrier of VjJ, indicating Fe^; 
has low mobility. 

FigureHKc) shows the structure of FeL;-Lip c . This anti- 
site pair has a formation energy of 0.51 eV. This value is 
independent of the choice of chemical potentials because 
the chemical potential term in the formation energy for- 
mula cancels out, cf. Eq. (0). Fe^-Li^ has a binding- 
energy of 0.44 eV; the distance between Fe^i and Lip e 
is 3.45 A, compared to 3.32 A between the lithium and 
iron sites. Finally, we find that 2Fe^--Vp r T has a forma- 
tion energy of 1.47—1.67 eV for reasonable choices of the 
atomic chemical potentials, and a binding energy of 1.25 
eV. With this high formation energy, the complex is un- 
likely to form in LiFePC>4, and is therefore not included 
in Fig. |31 Note that the formation energies of Fe^i-l^Li 
and 2Fe^ i -V^~ have the same dependence on the atomic 
chemical potentials (both contain the term — /ipc + 2/zlO 
and, hence, the same dependence on fio 2 ■ For any given 
set of chemical potentials, the formation energy of 2FeL ; - 
Vp~ is higher than that of Fc^-V^T by 1.11 eV. We also 
considered possible lithium and iron Frenkel pairs (i.e., 
interstitial- vacancy pairs), but these pairs are unstable 
toward recombination, probably because there is no en- 
ergy barrier or too small of a barrier between the vacancy 
and the interstitial. 

The above mentioned neutral defect complexes have 
also been studied by other research groups using cither 
interatomic-potential simulation :) 21 ' 22 or first-principles 
DFT calculations. 24 Islam et al. found that Fe^i-Lipe has 
a formation energy of 0.74 eV (or 1.13 eV if the two 



defects in the pair are considered as isolated defects) 
and a binding energy of 0.40 eV, and is energetically 
most favorable among possible native defects^! The re- 
ported formation energy is, however, higher than our cal- 
culated value by 0.23 eV. This difference may be due to 
the different methods used in the calculations. Fisher 
et al. reported a formation energy of 3.13 eV for Fe^i- 
V^,— which is much higher than our calculated value. 
Note, however, that Fisher et al. assumed the reaction 
FeO + 2Li£i-KFe£i + V u + Li 2 for the formation of Fe+- 
VjT which implies that LiFeP04 is in equilibrium with 
FeO and Li20. This scenario is unlikely to occur, as in- 
dicated in the Li-Fe-P-02 phase diagrams calculated by 
Ong et al.^i where equilibrium between these phases has 
never been observed. This may also be the reason that 
the formation energy of VjJJ reported by the same authors 
(4.41 eV)^ is much higher than our calculated values. 

Based on first-principles calculations, Malik et al. re- 
ported a formation energy of 0.515—0.550 eV for the anti- 
site pair Fe^i-Lipgri which is very close to our calculated 
value (0.51 eV). For 2Fel r Vp~ , the formation energy was 
reported to be of about 1.60—1.70 eV for no 2 ranging 
from —3.03 to —8.21 eV; 2 - 4 - which is also comparable to 
our results. Malik et al., however, obtained a much higher 
formation energy for Fe^-V^ , from about 3.60 to 5.10 eV 
for the same range of /zo 2 values^ This energy is much 
higher than that obtained in our calculations (0.36—0.56 
eV). Although we have no explanation for this discrep- 
ancy we observe that the calculated formation energies 
of 2Fe£ i -V F 2 ~ and Fe^-V^T in Malik et al.'s work have dis- 
tinct /io 2 -dependencies (see Figure SI in the Supporting 
Information of Refi^ 4 -), instead of having the same de- 
pendence on /zo 2 as we discussed above, indicating their 
scheme of accounting for the atomic chemical potentials 
differs from the standard procedure. 



V. TAILORING DEFECT CONCENTRATIONS 

It is important to note that the energy landscape pre- 
sented in Fig. [3] may change as one changes the atomic 
chemical potentials, i.e., synthesis conditions. The cal- 
culated formation energies are a function of four vari- 
ables /iLi, hfc, £*p, and /jo 2 , which in turn depend on 
each other and vary within the established constraints. 
A change in one variable leads to changes in the other 
three. In the following discussions, we focus on two 
"knobs" that can be used to experimentally tailor the 
formation energy and hence the concentration of different 
native defects in LiFeP04, and suppress or enhance cer- 
tain defects for targeted applications. One is /zo 2 , which 
can be controlled by controlling temperature and pres- 
sure and/or oxygen reducing agents. Lower ^o 2 values 
represent the so-called "more reducing environments," 
which are usually associated with higher temperatures 
and/or lower oxygen partial pressures and/or the pres- 
ence of oxygen reducing agents; whereas higher /zo 2 val- 
ues represent "less reducing environments." 41 The other 
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FIG. 5: Calculated formation energies of native point defects 
and defect complexes in LiFePCU, plotted as a function of 
Fermi level with respect to the VBM. The energies are ob- 
tained at /^o 2 — — 3.03 eV, and equilibrium with Fe203 and 
Fe7(P04)6 is assumed. 
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FIG. 6: Calculated formation energies of native point defects 
and defect complexes in LiFeP04, plotted as a function of 
Fermi level with respect to the VBM. The energies are ob- 
tained at fj,o 2 =—8.21 eV, and equilibrium with Fe2P and FesP 
is assumed. 



is the degree of lithium off-stoichiometry with respect to 
LiFeP04 exhibited through the tendency toward forma- 
tion of Li-containing or Fe-containing secondary phases 
in the synthesis of LiFeP04. As discussed previously, 
in the environments to which we refer as Li-excess (Li- 
deficient), the system is close to forming Li-containing 
(Fe-containing) secondary phases. 

Varying the Atomic Chemical Potentials. Let 
us assume, for example, Li-deficient environments and 
vary /io 2 from —3.03 (where LiFePCU first starts to 
form) to —8.25 eV (where it ceases to form)4i This 
amounts to choosing different cuts along the /^o 2 axis 
in Fig. [1] to give different two-dimensional polygons of 
LiFeP04 stability. Figure [5] shows the calculated forma- 
tion energies for fj,o 2 =— 3.03 eV, assuming equilibrium 
with Fe2 03 and Fe7(P04)6 (i.e., Li-deficient) which gives 
rise to (J>u=— 3.41, ^f c =— 3.35, and fJ,p=— 6.03 eV. Fig- 
ure [5] clearly shows changes in the energy landscape of 
the defects, compared to Fig. [31 The lowest energy point 
defects that determine the Fermi-level position are now 
p + and VjJ. Near ei?=0.59 eV where p + and VT, have 
equal formation energies, also has the lowest energy. 
This indicates that, under high /io 2 an d Li-deficient en- 
vironments, p + and are the dominant native point 
defects in LiFeP04 and are likely to exist in the form of 
the neutral complex V®- r Note that, with the chosen set of 
atomic chemical potentials, Lip c also has a low formation 
energy, very close to that of , indicating the presence 
of a relatively high concentration of Lip c . Similar to V^, 
Lip e can combine with p + to form Li Fe . However, since 
Li Fo has a higher formation energy and a smaller binding 
energy than V^; , only a small portion of Lip e is expected 
to be stable in form of Li Fc under thermal equilibrium 
conditions. Iron vacancies have the lowest energies in a 
wide range of the Fermi-level values as expected, given 



the very low iron chemical potential. 

Figure [5] shows the calculated formation energies for 
/io 2 =— 8.21 eV. The formation energies are obtained 
by assuming equilibrium with Fe2P and FesP (i.e., Li- 
deficient) which gives rise to //li=— 1-80, /iF C =— 0.24, and 
/ip=— 0.39 eV. We find that Fe^j and Lip e are now the 
dominant native point defects, pinning the Fermi level at 
£f=2.00 eV. The complex Fe^-Li^ has a binding energy 
of 0.44 eV, comparable to the formation energies of FeJ; 
and Li Fo (which are both 0.48 eV at e F =2.00 eV). This 
suggests that Fe^ and Lip e are likely to exist both in the 
form of FeLi-Lip G , but also as isolated point defects. With 
this set of atomic chemical potentials, we find that V^q^ 
has the lowest formation energy near the VBM, and Fe^ 
has a formation energy of 1.15 eV that, while very high, is 
lower than (1.93 eV) and Li Fo (1.98 eV), which were 
found to have lower formation energies under different 
conditions (c/. Figures 3 and 5). 

We also investigated the dependence of defect forma- 
tion energies on /XLi (and /iFc), i-e., Li-deficiency versus 
Li-excess, for a given /io 2 value. For ^q 2 —— 4.59 eV, for 
example, the results obtained at Points B and C in Fig.Q] 
show energy landscapes that are similar to that at Point 
A (Li-deficient), namely Fe^ and V" L 7 are the dominant 
native point defects in LiFeP04 and likely to exist in 
form of Fe^-V^- At Point D, where LiFeP04 is in equi- 
librium with L14P2O7 and Li3PC>4 (Li-excess), we find 
instead that Fe^; and Lip e are energetically most favor- 
able, and are likely to exist as FeLi-Lip . The calculated 
formation energy of p + is only slightly higher than that 
of Fe^i, indicating a coexisting high concentration of p + . 
The hole polarons in this case are expected to exist as 
isolated defects under thermal equilibrium since Li Fc has 
a relatively high formation energy (0.66 eV) and a small 
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binding energy (0.30 eV). Point E gives results that are 
similar to those at Point D. In contrast, when we choose 
Ho 2 =— 8.21 eV, we find that Fe^ and Lip c are the most 
energetically favorable defects regardless of the choice of 
phase-equilibrium conditions. 

Identifying the General Trends. We list in Table U 
the formation energies of the most relevant native point 
defects and defect complexes in LiFeP04, migration bar- 
riers of selected point defects, and binding energies of 
the defect complexes. The chemical potentials are cho- 
sen with representative /io 2 values and Li-deficient and 
Li-excess environments to reflect different experimental 
conditions. Specifically, these conditions represent equi- 
librium of LiFeP0 4 with (1) Fe 2 3 and Fe 7 (P0 4 ) 6 and 
(2) Li 3 Fe 2 (P0 4 )3 and Li 3 P0 4 , for /i O2 =-3.03 eV; (3) 
Fe 2 3 and Fe 3 (P0 4 ) 2 and (4) Li 4 P 2 7 and Li 3 P0 4 , for 
/io 2 =-3.89 eV; (5) Fe 2 3 and Fe 3 (P0 4 ) 2 (i.e., Point A 
in Fig. [Lj and (6) Li 4 P 2 7 and Li 3 P0 4 (i.e., Point D 
in Fig. ID, for /z 0a =-4.59 eV; (7) Fe 3 P and Fe 2 P and 
(8) Li 3 P0 4 and Fe 3 P, for /x 02 =-8.21 eV. Conditions 
(1), (3), (5), and (7) represent Li-deficient environments, 
whereas (2), (4), and (6) represent Li-excess. Under each 
condition, the formation energies for charged defects are 
taken at the Fermi-level position determined by relevant 
charged point defects: (1) e F =0.59 eV (where p + and 
have equal formation energies), (2) e F =0.58 eV (p + 
and Li"), (3) e F =0.79 eV (Fe+ and V^), (4) e F =0.66 
eV (p+ and Li Fc ), (5) e F =1.06 eV (Fe+ and Vjj), (6) 
e F =0.74 eV (Fe+ and Li Fc ), (7) e F =2.00 eV (Fe+ and 
Lip c ), and (8) e F =1.98 eV (Fe+ and Li Fc ). The results 
for other /zo 2 values are not included in Table Q] because 
they give results that are similar to those presented here. 
For example, the energy landscapes for /io 2 =— 3.25 eV 
and /zo 2 =— 7.59 eV are similar to those for ^o 2 =— 3.03 
eV and /Lto 2 =— 4.59 eV, respectively. 

In order to help capture the most general trends in 
the energy landscape of native defects in LiFeP0 4 in go- 
ing from high to low \xq 2 values and from Li-deficient to 
Li-excess environments, we plot in Fig. [7] the calculated 
formation energies of the most relevant native point de- 
fects and their lowest-energy complexes obtained under 
conditions (1) — (8). We find that, at a given Fermi-level 
position e F , the formation energy of Fe^; decreases as 
H0 2 decreases. This is because /ZFe increases more rapidly 
than /Lt Li does as /zo 2 decreases from —3.03 eV to —8.21 
eV. The formation energy of p + , on the other hand, is 
independent of the choice of atomic chemical potentials, 
and depends only on e F . At high fi 02 values, p + is lower 
in energy than Fe^, but then the two defects switch or- 
ders before fio 2 reaches —3.89 eV (under Li-deficient en- 
vironments) or —4.59 eV (Li-excess). The formation en- 
ergy of the dominant positive defects, p + and FeJ^, differs 
by as much as 1.3 eV for some sets of atomic chemical 
potentials. On the contrary, the dominant negative de- 
fects, V^T and Lipg, have comparable formation energies 
throughout the range of conditions. The largest forma- 
tion energy difference between the two defects is just 0.2 
eV. The formation energy of Lip is slightly lower than 



that of U L ~j" under Li-excess environments; whereas it is 
slightly higher under Li-deficient environments, except 
near /Lto 2 =— 8.21 eV where U L 7 is higher than Lip c . Both 
and Lip o have their formation energies increased as 
/U0 2 decreases. In going from high to low /j,o 2 values, the 
changes in the formation energy of Fe^; and that of U L ~[ 
and Lip e leads to a shift of the Fermi level from about 
0.6 eV above the VBM to 2.0 eV above the VBM. The 
variation in the calculated formation energy of p + as seen 
in Table U is a result of this shift. 

Under Li-deficient environments, we find that Fej^ and 
VjT are energetically most favorable over a wide range of 
fio 2 values, from —3.89 to —7.59 eV, and are likely to ex- 
ist in form of the neutral complex Fe^-VjJ . At the higher 
end in the range of fXo 2 values, from —3.03 to —3.25 eV, 
p + and U L 7 are the most favorable, and are likely to exist 
in form of the complex V®- r Finally, only at the lowest 
end of the fio 2 values, the most favorable point defects 
are Fe^; and Lip c , which may exist in form of the neu- 
tral complex Fe^-Lipe. Under Li-excess environments, 
Lip e dominates the negatively charged point defects in 
the whole range of /io 2 values. This makes p + and Lip e 
the most favorable point defects for fio 2 ranging from 
—3.03 to —3.89 eV, and FeLi an d Lip c the most favorable 
defects for fio 2 ranging from —4.59 to —8.21 eV. Note 
that, although the formation energy difference between 
VjJ and Lip e is small (less than 0.2 eV), the difference in 
their concentrations can still be significant, as indicated 
by the exponential dependence in Eq. ([1]). 

Overall, we find that the calculated formation ener- 
gies of the dominant native point defects are low, from 
about 0.3 to 0.5 eV for /io 2 from -3.03 to -8.21 eV 
(c/. Table HJ. With such low formation energies, the de- 
fects will easily form and occur with high concentrations. 
The dominant defects may be different, however, if one 
changes the experimental conditions during synthesis, as 
discussed above. This is consistent with the reported 
experimental data showing the presence of various na- 
tive defects in LiFeP0 4 samples*^— We note that there 
are several limitations inherent in our calculations. The 
first set of limitations comes from standard methodologi- 
cal uncertainties contained in the calculated formation 
enthalpies and phase diagrams as discussed in Ref. 41 . 
The second set comes from the calculation of defect for- 
mation energies using supercell models where supercell 
finite-size effects are expected^ Since applying approx- 
imations indiscriminately in an attempt to correct for 
finite-size effects tends to "overshoot" and makes the en- 
ergies even less accurate^ we did not include any cor- 
rections pertaining to such effects in our defect calcula- 
tions. A proper treatment of finite-size effects, however 
if applied, will lead to an increase in the calculated for- 
mation energy of the charged native point defects, and 
hence the binding energy of the neutral defect complexes. 
In spite of the limitations, the general trends discussed 
above should still hold true. And, our results therefore 
can serve as guidelines for tailoring the defect concentra- 
tions in LiFeP0 4 , and suppressing or enhancing certain 
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TABLE I: Calculated formation energies (E* ) and migration barriers (E m ) of the most relevant native point defects and defect 
complexes in LiFePOzi. (1) — (8) are the equilibrium conditions; see text. Binding energies (Eb) of the defect complexes (with 
respect to their isolated constituents) are given in the last column. The formation energy of 2Fe^ r Vp~ is high and thus the 
complex is not likely to form, but is also given here for comparison. 



Defect 


(1) 


(2) 


(3) 


E f (eV) 
(4) (5) 


(6) 


(7) 


(8) 


E m (eV) 


Constituents 


E b (eV) 


P + 


0.33 


0.32 


0.54 


0.40 


0.80 


0.49 


1.74 


1.72 


0.25 






Fe+ 


0.57 


0.63 


0.43 


0.55 


0.42 


0.48 


0.48 


0.48 








Vu 


0.33 


0.41 


0.43 


0.47 


0.42 


0.56 


0.53 


0.55 


0.32 






Lire 


0.39 


0.32 


0.52 


0.40 


0.53 


0.48 


0.48 


0.48 








vh 


0.32 


0.38 


0.62 


0.52 


0.88 


0.70 


1.93 


1.95 


0.32° 


v u + P + 


0.34, 0.42 6 


Li Pc 


0.42 


0.34 


0.76 


0.50 


1.03 


0.66 


1.92 


1.90 




Lipe + P + 


0.30 




0.41 


0.55 


0.37 


0.53 


0.36 


0.55 


0.52 


0.56 




Fe t + v u 


0.49 


Fe^-Lipe 


0.51 


0.51 


0.51 


0.51 


0.51 


0.51 


0.51 


0.51 




Fe^ ; + Li" 


0.44 


2Fe+-^ e - 


1.52 


1.66 


1.48 


1.64 


1.47 


1.66 


1.63 


1.67 




2Fe+ + Vl~ 


1.25 



a Lower bound, estimated by considering VS as a complex of V Li and p + and taking the higher of the migration energies of 
the constituents. 6 The value obtained in calculations using larger supercells (252 atoms/cell). 
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FIG. 7: Calculated formation energies of the low-energy positively and negatively charged point defects (i.e., p + , Fe^, V{^, 
and Lip e ) and their neutral defect complexes, plotted as a function of Fermi level with respect to the VBM, under different 
conditions: /j.o 2 — — 3.03, —3.89, —4.59, and —8.21 eV, and Li-deficient and Li-excess environments. Panels (a) — (h) correspond 
to conditions (1) — (8), see text. Only the complex that consists of the lowest-energy negatively and positively charged point 
defects are included. 



defects for targeted applications. 

VI. ELECTRONIC AND IONIC CONDUCTION 

Strictly speaking, lithium vacancies are only stable as 
V{^, and Vlj (which is, in fact, a complex of and p + ) 
cannot be considered the vacancies's neutral charge state. 
Likewise, lithium antisites, iron antisites, and iron vacan- 
cies also have one stable charge state only and occur as, 
respectively, Lip e , Fe^, and Vp~. Removing/ adding elec- 
trons from/to these stable point defects always results in 



defect complexes consisting of the point defects and small 
hole/electron polarons, as presented in the previous sec- 
tions. The fact that small polarons can be stabilized, 
both in the absence and in the presence of other native 
defects is necessarily related to the electronic structure 
of LiFeP04 where the VBM and CBM consist predom- 
inantly of the highly localized Fe 3d states. Combined 
with the fact that charged native defects have negative 
formation energies near the VBM and CBM (c/. Figures 
3, 5, and 6), our results therefore indicate that native de- 
fects in LiFePCU cannot act as sources of band-like hole 
and electron conductivities. These defects will, however, 
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act as compensating centers in donor-like doping (for the 
negatively charged defects) or acceptor-like doping (the 
positively charged defects). The electronic conduction in 
LiFePC>4 thus occurs via hopping of small hole polarons. 
This mechanism, in fact, has been proposed for LiFePC>4 
in several previous works.— ~— Zaghib et al^ found ex- 
perimental evidence of intra-atomic Fe 2+ — Fe 3+ transi- 
tions in the optical spectrum of LiFeP04, which indi- 
rectly confirms the formation of small hole polarons. The 
activation energy for electronic conductivity in LiFeP04 
was estimated to be 0.65 eV^ comparable to that of 
0.55—0.78 eV reported by several other experimental 
groups ^ 12 ' 51 

In order to compare our results with the measured ac- 
tivation energies, let us assume two scenarios for hole 
polaron hopping in LiFeP04. In the first scenario, we 
assume self-diffusion of free p + defects. The activation 
energy E a for this process is calculated as the summa- 
tion of the formation energy and migration barrier of p + , 
where the former is associated with the intrinsic concen- 
tration and the latter with the mobility, 52 

E a . = Ef(p+)+E m (p+), (7) 

which gives i? a =0.57 eV, if E*(p + ) is taken under the 
most favorable condition where p + has the lowest for- 
mation energy (cf. Table [I]). In the second scenario, we 
assume that p + and are formed via the formation of 
the neutral complex V^, similar to cases where defects 
are created via a Frenkel or Schottky mechanism^ At 
high temperatures, the activation energy for the diffusion 
of p + can be calculated as 

E a = ±Ef(V° i ) + E m (p+), (8) 

which results in E a =0Al eV, assuming the condition 
where V® { has the lowest formation energy (cf. Ta- 
ble |T|. The lower bound of the activation energy for 
polaron conductivity is therefore 0.41—0.57 eV. This 
range of E a values is comparable to that obtained in 
experiments 4 n 12 i 13 i 51 

Among the native defects, is the most plausible 
candidate for ionic conduction in LiFeP04, because of its 
low formation energy and high mobility. Using formulae 
similar to Eqs. (7) and (8), we estimate the lower bound 
of the activation energy for self-diffusion of VjJ along the 
&-axis to be 0.48—0.65 eV. The reported experimental 
value is 0.54 eV— or 0.62—0.74 eV, 51 obtained from ionic 
conductivity measurements carried out on LiFeP04 sin- 
gle crystals, which is comparable to our calculated value. 
The diffusion of V^, however, may be impeded by other 
native defects or extrinsic impurities that have lower mo- 
bility. The presence of Fe^ in LiFeP04 has been reported 
and the defect is believed to reduce the electrochemi- 
cal performance of the material by blocking the lithium 
channels*^— Indeed, our results also show that Fe^ oc- 
curs with a high concentration under various conditions 
(cf. Table IJ) and has low mobility. Whether Fe^; is sta- 
ble in form of the neutral complex Fe^-V^ or Fe^j-Lipg 



as suggested in several experimental works ; 15 ! — how- 
ever, depends on the specific conditions under which the 
samples are produced. 

What is most fascinating about our results is that one 
can suppress Fej^ by adjusting suitable experimental con- 
ditions during synthesis. In fact, p + -rich and FeLj-free 
LiFeP04 samples, which are believed to be desirable for 
high intrinsic electronic and ionic conductivities, can be 
produced if one maintains high /io 2 values (cf. Fig. [7]). 
Of course, /j,o 2 should not be so high that LiFePC>4 be- 
comes unstable toward forming secondary phases. Al- 
though LiFePC>4 cannot be doped p-type or n-type as 
discussed earlier, the incorporation of suitable electri- 
cally active impurities in the material can enhance the 
electronic (ionic) conductivity via increasing the concen- 
tration of p + (U L T). These impurities, if present in the 
samples with a concentration higher than that of the 
charged native defects, can shift the Fermi level^r— and 
hence lower the formation energy of either p + or VjJ . A 
decrease in the formation energy of p + , however, may 
result in an increase in the formation energy of and 
vice versa. For example, impurities with positive effec- 
tive charges (i.e., donor- like doping) may shift the Fermi 
level to the right (cf. Fig. [7]), resulting in an increased (de- 
creased) formation energy of p + (VjJ). Impurities with 
negative effective charges (i.e., acceptor-like doping), on 
the other hand, may produce the opposite effects, namely, 
decreasing (increasing) the formation energy of p + (V^). 
An enhancement in both electronic and ionic conductiv- 
ities would, therefore, require a delicate combination of 
defect-controlled synthesis, doping with suitable electri- 
cally active impurities, and post-synthesis treatments. 
An example of the latter would be thermal treatment 
which, in fact, has been reported to cause lithium loss in 
LiFePC>4 and lower the activation energy of the electrical 
conductivity^ 



VII. CONCLUSION 

In summary, we have carried out comprehensive first- 
principles studies of native point defects and defect com- 
plexes in LiFePC>4. We find that lithium vacancies, 
lithium antisites, iron antisites, and iron vacancies each 
have one stable charge state only and occur as, respec- 
tively, V^, Lip e , Fe^i, and . The removal/addition of 
electrons from/to these stable native point defects does 
not result in a transition to other charge states of the 
same defects, but instead generates small hole/electron 
polarons. The fact that small polarons can be stabilized, 
both in the presence and in the absence of other native 
defects, is necessarily related to the electronic structure 
of LiFcPC>4. Our analysis thus indicates that native de- 
fects in the material cannot act as sources of band-like 
electron and hole conductivities, and the electronic con- 
duction, in fact, proceeds via hopping of small hole po- 
larons (p + ). The ionic conduction, on the other hand, 
occurs via diffusion of lithium vacancies. 
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Among all possible native defects, p + , V^, Lip c , and 
Fe^; are found to have low formation energies and are 
hence expected to occur in LiFePC>4 with high concen- 
trations. The dominant point defects in the samples 
are likely to exist in forms of a neutral defect complex 
such as Ligfc, Fe^-V^, or Fe^-Li Fc . The energy 
landscape of these defects is, however, sensitive to the 
choice of atomic chemical potentials which represent ex- 
perimental conditions during synthesis. This explains 
the conflicting experimental data on defect formation in 
LiFePOzt. Our results also raise the necessity of having 
prior knowledge of the native defects in LiFeP04 sam- 
ples before any useful interpretations of the measured 
transport data can be made. We suggest that one can 
suppress or enhance certain native defects in LiFeP04 
via tuning the experimental conditions during synthesis, 
and thereby produce samples with tailored defect concen- 



trations for optimal performance. The electrical conduc- 
tivity may be enhanced through increase of hole polaron 
and lithium vacancy concentrations via a combination of 
defect-controlled synthesis, incorporation of suitable elec- 
trically active impurities that can shift the Fermi level, 
and post-synthesis treatments. 
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